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Processes mediated by two virtual low-energy photons contribute quite significantly to the energy of hydro- 
genic 5 states. The corresponding level shift is of the order of (a/ir) 2 (Za) 6 m e c 2 and may be ascribed to a 
two-loop generalization of the Bethe logarithm. For 15 and 25 states, the correction has recently been eval- 
uated by Pachucki and Jentschura [Phys. Rev. Lett. 91, 113005 (2003)]. Here, we generalize the approach to 
higher excited 5 states, which in contrast to the 15 and 25 states can decay to P states via the electric-dipole 
(El) channel. The more complex structure of the excited-state wave functions and the necessity to subtract P- 
state poles lead to additional calculational problems. In addition to the calculation of the excited-state two-loop 
energy shift, we investigate the ambiguity in the energy level definition due to squared decay rates. 
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I. INTRODUCTION 

Both the experimental and the theoretical study of radia- 
tive corrections to bound-state energies have been the subject 
of a continued endeavor over the last decades (for topical re- 
views see fjl 0, IH 01). Simple atomic systems like atomic 
hydrogen, and heliumlike or lithiumlike systems, provide a 
testbed for our understanding of the fundamental interactions 
of light and matter, including the intricacies of the renormal- 
ization procedure and the complexities of the bound-state for- 
malism. One of the historically most problematic corrections 
for bound states in hydrogenlike systems is the two-loop self- 
energy (2LSE) effect (relevant Feynman diagrams in Fig.^, 
and this correction will be the subject of the current paper. 
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FIG. 1 : Two-photon processes may be interpreted in terms of Feyn- 
man diagrams. The double line denotes the bound-electron propaga- 
tor. In the figure, we display the crossed-loop (A), the rainbow (B), 
and the loop-after-loop (C) diagram. 

Regarding self-energy calculations, two different ap- 
proaches have been developed for hydrogenlike systems: (i) 
the semianalytic approach, which is the so-called Za ex- 
pansion and in which the radiative corrections are expressed 
as a semianalytic series expansion in the quantities Za and 



ln[(Za) -2 ], and (ii) the numerical approach, which avoids 
this expansion and leads to excellent accuracy for systems 
with a high nuclear charge number. Over the last couple of 
years, a number of calculations have been reported that profit 
from recently developed numerical algorithms and an im- 
proved physical understanding of the problem at hand. These 
have led to numerical results even at low nuclear charge num- 
ber 
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FIG. 2: (color online). Hydrogenic 5 levels have the same (radial) 
symmetry properties as the ground state. The wave function ip(r) 
of an 5 level therefore depends only on the radial coordinate r = 
\Jx 2 + y 2 + z 2 . The function fes(x,y) = J^dz \ipes(x, y, z)\ 2 
is positive definite and constitutes effectively an integrated projec- 
tion of the 65 electron probability density onto the x-y plane. In- 
deed, we plot here the natural logarithm of this function, which 
is ln[/6s(x, y)\, as a function of x £ [— 40 aBohr, 40 aeohr] and 
y 6 [—40 aeohr, 40 aeohr]- Here, aBohr denotes the Bohr radius 
a B ohr = h/(amc) = 0.529 177 2108(18) x 10" 10 m @|. 

Within approach ( i), a number of calculations have recently 
been reported which rely on a separation of the energy scale(s) 
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of the virtual photon(s) into high- and low-energy domains 
(see, e.g., 0, Chap. 123], O Chap. 7] and (H). This has 
recently been generalized to two-loop corrections lll2l ll3Tl . 
Also, there have been attempts to enhance our understanding 
of logarithmic corrections in higher orders of the Za expan- 
sion by renormalization-group techniques fl4lfl5ll . 

In the current paper, we discuss the evaluation of a specific 
two-loop correction, which can quite naturally be referred to 
as the two-loop generalization of the Bethe logarithm, for 
higher excited S states (see also Fig. 0. The calculation is 
carried out for the dominant nonlogarithmic contribution to 
the two-loop self-energy shift of order a 2 (Za) 6 m c 2 , where 
771 c 2 is the rest energy of the electron (m is the rest mass), 
a is the fine-structure constant, and Z is the nuclear charge 
number. 

The two-loop Bethe logarithm is formally of the same or- 
der of magnitude [a 2 (Za) 6 mc 2 ] as a specific set of two- 
loop corrections which are mediated by squared decay rates 
and whose physical interpretation has been shown to be lim- 
ited by the predictive power of the Gell-Mann-Low theo- 
rem on which bound-state calculations are usually based 1 1 6 1 . 
For excited nS states (n > 3), the problematic squared 
decay rates lead to an ambiguity which we assign to the 
two-loop Bethe logarithm as a further theoretical source of 
uncertainty. Beyond the order of a 2 (Za) 6 mc 2 , the def- 
inition of an atomic energy level becomes ambiguous, and 
the evaluation of radiative corrections to energy levels has 
to be aug mented by a more complete theory of the line 
shape (11 E3 EE Sim IH, with the IS state being the 
only true asymptotic state and therefore — in a strict sense — 
the only valid in- and out-state in the calculation of 5-matrix 
type amplitudes IT3l . Further interesting thoughts on ques- 
tions related to line shape profiles can be found in 11231 12411 . It 
had also been pointed out in Sec. VI of ll25ll that the asymme- 
try of the natural line shape has to be considered at the order 
a 8 . 

It is tempting to ask how one may intuitively understand 
the slow convergence of the Za expansion of the two-loop 
energy shift. As pointed out in y|, terms of different order in 
the Za expansion have rather distinct physical origins. In the 
order a 2 (Za) 4 mc 2 , there are two corrections, both arising 
from hard (high-energy) virtual photons. These correspond, 
respectively, to the infrared convergent slope of the two-loop 
electron Dirac form factor and to a two-loop anomalous mag- 
netic moment correction. The term of order a 2 (Za) 5 mc 2 
may also be computed without any consideration of low- 
energy virtual quanta yl [Tl], |26l [271 12811 . The terms of order 
a 2 (Za) 6 771 c 2 are not the leading terms arising from the two- 
loop self-energy shifts. Logarithmic correction terms of order 
a 2 (Zafh^Za)- 2 ]mc 2 (i = 1, 2, 3) have been consid- 
ered in (illl. It is only at the order of a 2 (Za) 6 m c 2 that 
the low-energy virtual photons begin to contribute to the hy- 
drogenic energy shift(s) of S states. They do so quite signif- 
icantly, enhanced by the triple logarithm (i = 3) and a sur- 
prisingly large coefficient of the single logarithm (i = 1), as 
shown in 111211 . It is therefore evident that we need to gain 
an understanding of all logarithmic and nonlogarithmic terms 
(i = 0,1,2,3) of order a 2 (Za) 6 \n[(Za)- 2 }mc 2 before 



any reliable prediction for the two-loop bound-state correc- 
tion can be made even at low nuclear charge number. 

An interesting observation can be made based on the fact 
that the imaginary part of the (nonrelativistic) one-loop self- 
energy gives the leading-order contribution to the El one- 
photon decay width of excited states l30tl . Analogously, it 
is precisely the imaginary part of the nonrelativistic two-loop 
self-energy which corresponds to the two-photon decay rate 
of the 2S state. The 25 two-photon decay rate is of the order 
of a 2 (Za) 6 m c 2 (see, e.g., JHHHIHIsl). From a nonrel- 
ativistic point of view, the scaling a 2 (Za) 6 m c 2 can be seen 
as some kind of "natural" order for the two-loop effect. It is 
the first order in which logarithms of Za appear and the first 
order in which a matching of low- and high-energy contribu- 
tions is required. This is also reflected in the properties of the 
two-photon decay. 

This article is organized as follows. In Sec. |ll] we review 
the status of known two-loop self-energy corrections. The for- 
mulation of the problem in nonrelativistic quantum electro- 
dynamics (NRQED) and calculation is discussed in Sec. 11111 
Squared decay rates are the subject of Sec. II VI and further 
contributions to the self-energy in the order of a 2 (Za) 6 m c 2 
are discussed in Sec. [V] Finally, conclusions are drawn in 
Sec.lVll 

II. KNOWN TWO-LOOP SELF-ENERGY COEFFICIENTS 

We work in natural units (h = c = £q = 1), as is cus- 
tomary in QED bound-state calculations. The (real part of 
the) level shift of a hydrogenic state due to the two-loop self- 
energy reads 

For the two-loop self-energy (2LSE) diagrams (see Fig.Q, the 
first terms of the semianalytic expansion in powers of Za and 
ln[(Za)- 2 ] read 

H(Za) = B^ SE) + (Za) B^ SE) 

+ (Za) 2 {Bg LSE) \n\Za)- 2 + B^ SE) \n 2 (Za)- 2 

+B { 6 f SE) \n(Za)~ 2 + B^ LSE) } . (2) 

The function H(Za) is dimensionless. We ignore unknown 
higher-order terms in the Za expansion and focus on a spe- 
cific numerically large contribution to B^ SE ^ given by the 
two-loop Bethe logarithm. We also keep the upper index 
(2LSE) in order to distinguish the two-loop self-energy con- 
tributions to the analytic coefficients from the self-energy 
vacuum-polarization (SEVP) effects il2ll35ll and the vacuum- 
polarization insertion into the virtual photon line in the one- 
loop self-energy (SVPE). By contrast, the sum of these effects 
carries no upper index, according to a convention adopted pre- 
viously in IfZ l35ll . It has been mentioned earlier that B4Q 
and B50 are purely relativistic effects mediated by hard vir- 
tual photons. The coefficient B40 in Eq. (|2ji involves a Dirac 
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and a Pauli form factor correction and reads 13 



B 



(2LSE) 
40 



(nS) = ~ - jjjj C(2) + 9 ln(2) C(2) - \ C(3) ; 

(3) 

the numerical value is 1.409 244. The first relativistic correc- 



tion B, 



(2LSE) 
50 



nS) is known to have a rather large value |2 



B^f SE \nS) = -24.2668(31). 
The triple logarithm in the sixth order of Za reads, 

B 63 (nS)=B^ SE \nS) = -^. 



(4) 



(5) 



It has meanwhile been clarified |37| |3J, |3JJ |40J] that the total 
value of this coefficient is the result of subtle cancellations 
among the different diagrams displayed in Fig.^ The double 
logarithm for nS is given by 



B, 



B. 



(2LSE) 
62 

(2LSE) 



(15) = S " 4 M2 



27 



-0.639 669, (6) 



(nS) = B 



(2LSE) 
62 



(IS) 



16 



T + t^t - - - ln(n)+*(n) + C 
4 An 2 n w w 



(7) 



where ^ denotes the logarithmic derivative of the gamma 
function, and C = 0.577216 ... is Euler's constant. 

The result for Ski, restricted to the two-loop diagrams in 
Fig.E reads (HQ 



b£ lse) (IS) 



5 221 875 9 , „ 

— C(2) + - C(2) In 2 



1296 
9 



C(3) 



152 

~2T 



In 2 



40, ,„ 
— In 2 2 
9 



+ -N(1S) 
49.838 317. 



(8) 



B, 



(2LSE) 
61 



(nS) 



Bi 2LSE \lS) + -[N(nS)-N(lS)] 



80 32 , „ 
27-¥ ln2 



1 

An 2 



- - ln(n) + *(n) + C 
n 



(9) 



We correct here a calculational error in Eq. (7a) of 
Ref. 1 3 511 where a result of 49.731651 had been given for 
B^ SE) (1S). However, even with this correction, the re- 
sult for B<; 2LSE) (15) given in Eq. is incomplete because 
it lacks contributions from two-Coulomb-vertex diagrams. 
These diagrams give rise to an effective interaction propor- 
tional to E 2 in the NRQED Hamiltonian and will be dis- 
cussed in l4lll . The additional contribution to B^ SE \lS) 
does not affect the n-dependence of this coefficient as indi- 
cated in Eq. l|9}, nor does it affect the calculation of the two- 
loop Bethe logarithms presented in this article. 



Numerical values of N(nS) are given in Eq. (12)] for 
= 1,...,8: 



N(1S) 


= 17.855 672(1), 


(10a) 


N(2S) 


= 12.032 209(1), 


(10b) 


N(3S) 


= 10.449 810(1), 


(10c) 


N(4S) 


= 9.722413(1), 


(10d) 


N(5S) 


= 9.304114(1), 


(10e) 


N(6S) 


= 9.031832(1), 


(lOf) 


N(7S) 


= 8.840123(1), 


(10g) 


N(8S) 


= 8.697639(1). 


(lOh) 



III. TWO-LOOP PROBLEM IN NRQED 

Historically, one of the first two-photon problems to be 
tackled theoretically in atomic physics has been the two- 
photon decay of the metastable 25 level which was treated 
in J43II44I1 . It is this decay channel which limits the lifetime 
of the 25 hydrogenic state. We have ll3lll 



= r = 



.229 Z 6 s~ 1 



= 1.310 Z 6 Hz. 



(11) 



The numerical prefactors of the width is different when ex- 
pressed in inverse seconds and alternatively in Hz. The fol- 
lowing remarks are meant to clarify this situation as well as 
the entries in Table [H] below. In order to obtain the width 
in Hz, one should interpret the imaginary part of the self- 
energy 1 30 1 as T/2, and do the same conversion as for the 
real part of the energy, i.e. divide by h, not fi. This gives the 
width in Hz. The unit Hz corresponds to cycles per second. 

In order to obtain the lifetime in inverse seconds, which is 
radians per second, one has to multiply the previous result by 
a factor of 2ir, a result which may alternatively be obtained 
by dividing T — i.e. the imaginary part of the energy, by h, not 
h. The general paradigm is that in order to evaluate an energy 
in units of Hz, one should use the relation E = hv, whereas 
for a conversion of an imaginary part of an energy to the in- 
verse lifetime, one should use T = hr^ 1 . As calculated in 
Refs. i3lll32il . the width of the metastable 25 state in atomic 
hydrogenlike systems is 8.229 Z 6 s _1 (inverse seconds). At 
Z = 1, this is equivalent to the "famous" value of 1.3 Hz 
which is nowadays most frequently quoted in the literature. 
The lifetime of a hydrogenic 25 level is thus 0.1215 Z~ 6 s. 
This latter fact has been verified experimentally for ionized 
helium 

We now briefly recall the expression for the two-photon de- 
cay involving two emitted photons with polarization vectors 
£1 and £2, in a two-photon transition from an initial state \<f>i) 
to a final state \4>f). The two-photon decay width T is given 
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by [see, for example, Eq. (3) of Ref. 13 111 1 

^max 

7T J 



A_ a 2 
27 



1 



H - E + u} 2 

e I 

H - E + uj 1 



where uj 2 = w max — io\ and w max = E — E' is the maximum 
energy that any of the two photons may have. The Einstein 
summation convention is used throughout this article. Note 
the identity I4IH49I1 



V 



1 



V 



m H — E + uji m 



f 



V 



1 



P' 



m H — E + u>2 m 



-u>i lo 2 m 



H - E + wi 



H -E + u} 2 



which is valid at exact resonance lo\ + uj 2 = Ei — Ej. 
This identity permits a reformulation of the problem in the 
velocity-gauge as opposed to the length-gauge form. 



(12) I n a number of cases, the formulation of a quantum elec- 
trodynamic bound-state problem may be simplified drasti- 
cally when employing the concepts of an effective low-energy 
field theory known as nonrelativistic quantum electrodynam- 
ics l5(ill . The basic idea consists in a correspondence be- 
tween fully relativistic quantum electrodynamics and effective 
low-energy couplings between the electron and radiation field, 
which still lead to ultraviolet divergent expressions. However, 
the ultraviolet divergences may be matched against effective 
high-energy operators, which leads to a cancellation of the 
cut-off parameters. Within the context of the one-loop self- 
energy pro bleir L a specialized approach has been discussed 
in fill l25l I5H l52ll . The formulation of the fwo-loop self- 
energy problem within the context of nonrelativistic quantum 
electrodynamics (NRQED) has been discussed in 11211 . We 
denote by the Cartesian components of the momentum op- 
erator p = — i V. The expression for the two-loop self-energy 

(13) shift reads E3 El 



A^nrqed = - ( Q 2 \ ) / dwiwi ( duj 2 uj 2 lip 1 — \r, P 3 — — ■ P l — ^- P> 

\ 3 7r m z J J J l\ H - E + uji H - E + loi+uj 2 H - E + uj 2 t 

1 1 1 

p p J p J 



2 \ H - E + oji H - E + oj 1 +oj 2 H -E + u>i 



— \p" — — p 3 — — p J 



2 \ H - E + uj 2 H - E + wx + uj 2 H - E + uj 2 

i—p-f ' V^_! A 



H-E + uji 1 \H-EJ 1 H-E + uj 2 ' 

1 A / • / 1 \ 2 \ 1 / 1 A / ■ / 1 



! \ pl H-E + ./ j \ p3 { h-E + . 2 ) P 7-2\^ H-E + u/ l \f [ H-E + .J P " ] 

* — — A - U — I — A m U — - — A) .(i4) 

H-E + ui H-E + U2 / u>i + lu 2 \ H — E + lu 2 / ui + uj 2 \ H - E + u x ' J 



All of the matrix elements are evaluated on the reference state 
\4>), for which the nonrelativistic Schrodinger wave function 
is employed. The Schrodinger Hamiltonian is denoted by H, 
and E = —(Za) 2 m/(2n 2 ) is the Schrodinger energy of the 
reference state (n is the principal quantum number). 

Expressions H2\ and dl 31 now follow in a natural way as 
the imaginary part generated by the sum of the first three terms 



in curly brackets in Eq. J14t . Specifically, the poles are gener- 
ated upon u>2 -integration by the propagator 



H — E + lo\ + lo 2 ' E; — E + wi + uj 2 

3 J 

when E — Ej = u>\ + ui 2 , which is just the energy conserva- 
tion condition for two-photon decay. Of course, other terms 
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in Eq. dl4> . not just the first three in curly brackets, may also 
generate imaginary parts (especially if the reference state is 
an excited state, and one-photon decay is possible). The cor- 
responding pole terms must be dealt with in a principal-value 
prescription if we are interested only in the real part of the 
energy shift. For P states and higher excited S states, the 
remaining imaginary parts find a natural interpretation as ra- 
diative correction to the one -photon decay width 1 54] . 
From here on we scale the photon frequencies cj± t 2 by 



<4 



UJf. 



(Za) 2 m 



1,2. 



(16a) 



This scaling, which is convenient for our numerical calcula- 
tions, (almost) corresponds to a transition to atomic units [but 
with (Za) 2 m = 1 instead of a unit Rydberg constant]. The 
momentum operator is scaled as 



P 



P 



Za m 



(16b) 



and becomes a dimensionless quantity. The Schrodinger 
Hamiltonian is scaled as 



H -> H' 



H 



(16c) 



(Za) 2 m 

The binding energy of the reference state receives a scaling as 



E —> E' 



E 



(Za) 



1 

Zn 2 



and is from now on also a dimensionless quantity (n is the 
principal quantum number). The scaled, dimensionless radial 



J 



coordinate is obtained as 

r — > i 

The scaled Green function 

GV) = 



Za to r . 



1 



E' — H' — oj> 
is also dimensionless. Finally, the quantity 



(16e) 



(16f) 



j) 01 



E> - E'. 



(16g) 



is the reduced Green function where the reference state \<f>) 
is excluded from the sum over intermediate states. The 
(scaled dimensionless) Schrodinger Hamiltonian is then given 
as H' =p' 2 /2— 1/r'. Scaled quantities will be used from here 
on until the end of the current Section|^ and we will denote 
the scaled, dimensionless quantities by primes, for absolute 
clarity of notation. (Note that in Ref. [55||, the correspond- 
ing scaled quantities were denoted by the same symbol as the 
dimensionful quantities.) As indicated in l55l Eq. (5)], the ex- 
pression dl4t can be rewritten in terms of the scaled quantities 
as 



(16d) A^nrqed = g (^) (Zafm j dcj[ J d^/K,^), 



(17) 

where the (dimensionless) function f(u)[,u> 2 ) is defined as 
[see also Eq. J14H 



(p H G'(u[)p' j G'fai + cj 2 )p H G'(J 2 )p'i) + - (p H G'^p'i G'(uj[ + J 2 )p'^ G'(uj[)p" 1 ) 



+ - ( P H G'(uj' 2 )p' j G'(cj[ + J 2 )p'i G>(J 2 )p*) + (p H G'( U [)p" G'U0)P' 3 G'(u) 2 )p H ) 



+ {p" l G'(u' 1 )G'(u' 2 )p' 1 ) 



1 



1 



- </G'(^)/) - -r—- (p H G'(uj[)p' 1 ) 



(18) 



r 



In S, the corresponding Equation (5) has a typographical 
error: the term to (p n G'(wJ) G'(cu 2 )p n ') should have a plus 
instead of a minus sign (seventh term in the square brack- 
ets). In particular, the scaling d!6i leads to a disappear- 
ance of the powers of Za when considering the expression 
/ duj i u[ J duj' 2 J 2 f(oj[,uj 2 ). 

First, we fix lo' x and integrate over uj' 2 . The subtraction pro- 
cedure is as follows. We need to subtract the contribution 
from the following terms that lead to divergent expressions 



as uin 



oo. We therefore expand f(w'\,io 2 ) for large cu' 2 at 



fixed uj' 1 . The asymptotics read |55| 



6K) 



(19) 



where the further terms in the expansion of f(ui[,ui' 2 ) for 
uj 2 — ► oo lead to convergent expressions when integrated over 



6 



uj' 2 in the region of large u) 2 - The leading coefficient is 



(H' -E' + uj' 1 f 



and the second reads 



(20) 



(21) 



where by Sw we denote the first-order correction to the quan- 
tity in curly brackets obtained via the action of the scaled, 
dimensionless, local potential 



W 



= 7rd {3 \r') 



(Za) 3 m 3 

i.e., by the replacements [see Eq. (I16A 1. 

H' -> H' + W, 

|0) -> |0> + |<ty>, 
E 1 -> E' + SE'. 

Here 



(22) 

(23a) 
(23b) 
(23c) 

(24) 



The correction J2 1> has been calculated for excited S states in 
Ref. IH. 



numerical data 

large-u/, asymptotics 




FIG. 3: (color online). Plot of the large-wj asymptotics of g [see 
Eq. <29H against numerical data obtained for g3s{w[) in the range 
u>[ £ [20, 180]. The numerical data are scaled by a factor of 3 3 = 
27. See also Tableland Eq. <25a> where g is defined. Dimensionless 
quantities are displayed in the figure [this statement relates to both 
the abscissa as well as the ordinate axis, see alsoEq. fifll. 

We are interested in evaluating the constant term g{u>'±) in 
the integral of f{u>i,uj' 2 ) in the range lu 2 G (0, A) at fixed uj^ 
for large A: 

/ dc^/ViM) = a(u}[) A + 6(^i) \nA + g(Lu[) , (25a) 



where we neglect terms that vanish as A — > oo. This equation 
provides an implicit definition of g(co[) as the constant term 
which results in the limit A — > 0. The constant term may be 
evaluated as 



where 



li = 



li = 



M 



d^2 J '(u'i,u' 2 ) , 



du)' 



M 



/(«>i,<4) - a(aj[) - 



1 3 = -a(u[)M- &K) InM, 



(25b) 



(26a) 



, (26b) 



(26c) 



with arbitrary M [the result for g(cj[) is independent of M]. 
Sample values of the g-function for nS states are given in Ta- 
ble|J The sign of the23-term (cf. I55I Eq. (8c)]) is determined 
by the necessity of subtracting the integral of the subtraction 
term [second term in the integrand of Eq. J26bH . at the lower 
limit of integration M. In both the ui[ as well as the oj' 2 in- 
tegrations, there is a further complication due to bound-state 
poles (P states) which need to be considered for higher ex- 
cited nS states (n > 3). In the current section, we completely 
ignore the imaginary parts and carry out all integrations with a 
principal-value prescription. Idem est, we use the prescription 
(M > a) 



M 



da/ 



In 



M 



(27) 



For double poles, which originate from some of the terms in 
Eq. H4\ . the appropriate integration prescription is as follows: 



M 



dcj' 



1 



M 



{u> - a) 2 a (a - M) 



(28) 



Even if M > a, this prescription leads to a finite result which 
is real rather than complex. The same result can also be ob- 
tained under a symmetric deformation of the integration con- 
tour into the complex plane. Analogous integration prescrip- 
tions have been used in I25I IBUl . Double poles normally lead 
to nonintegrable singularities and give rise to serious concern. 
It is therefore necessary to ask how these terms originate in 
the context of the current calculation. To answer this ques- 
tion it is useful to remember that expression dl4> is obtained 
by perturbation theory in powers of the nonrelativistic QED 
interaction Lagrangian; an expansion in powers of the inter- 
action is, however, not allowed when we are working close 
to a resonance of the unperturbed atomic Green function — 
i.e., close to a bound-state pole. The double poles incurred by 
this expansion find a natural a posteriori treatment by the pre- 
scriptions J27t and d28l > above. In general, double poles as en- 
countered here and previously in 1251 15111 originate whenever 
we work with (i) excited states which can decay through El 
one-photon emission and (ii) propagators are perturbatively 
expanded near bound-state poles. 
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The two-loop Bethe logarithm, which is equal to the low- 
energy contribution B 1 ^ (nS) to the coefficient -Bgo L [see 
Eq. ©], can be obtained by considering 



(31) 



and subtracting all terms that diverge as A — > oo, as given by 
the leading asymptotics in Eq. ( I29> . Specifically, the integra- 
tion procedure is as follows. We define the two-loop Bethe 
logarithm as 



b L (nS) =n 3 (J1+J2 + J3) 



(32) 



FIG. 4: (color online). Large-tvJ asymptotics of g plotted against 
numerical data for the 65* state in the range uj[ £ [20, 180]. The nu- 
merical data are scaled by a factor of 6 3 = 216. Explicit numerical 
sample values for ges(u'i) can also be found in Table|l| The appar- 
ent similarity of Figs.[3]and|4]is reflected in the scaled entries in this 
Table. The difference between the numerical data (solid line) and the 
asymptotics (dashed line) is negative. The formula for the large-c^ 
asymptotics of g is given in Eqs. 1291 and J30> . The difference be- 
tween the numerical data and the asymptotics gives rise to a negative 
contribution to the integral J2 defined according to Eq. I33bi and to 
a negative value for the Beo -coefficient [see Eq. 1521 below]. The 
scaled, primed quantities plotted along both the abscissa as well as 
the ordinate axis are dimensionless [see also Eq. 1161 1. 



where 



Ji = 

J2 = 



N 



du>[ 



N 



+C 



InK) 



D 



ln 2 K) 



F 



(33a) 



7t A \nu)[ 

n 6 



D 



(33b) 



The leading terms in the asymptotics of g(to[) for large u[ 
read 



where 



1 



9(^i) = Q 



-E 



A \tluj[ + B + C 
ln 2 K) p 1 



InK) 



D 



A = -4, 

B = 2 [In 2 - 1 - \nk (nS)] , 

C = 4V2, 

D = 4\/2 (2 (la 2-1) -tt) , 

£7 = 1, 



F = 8 + - N{nS) + 5tt 2 



1 



(29) 

(30a) 
(30b) 
(30c) 
(30d) 
(30e) 

(30f) 



The higher-order terms in the large-t^ expansion, which are 
ignored in Eq. ( I29> . lead to convergent expressions in the 
problematic integration region ui[ — > 00. Explicit numerical 
values for N(nS) are given in Eq. J 1 0b . For 3S and 65 states, 
numerical data for g are compared to the leading asymptotics 
in Figs.[5]and|4] 



J 3 = AN(lnN- 1) + BN 

+C VN (In N - 2) + 2D Vn 
+2EVN (8 + (In TV - 4) ln7V) 
+F\nN . 



Again, in analogy to the integration prescription in Eqs. J25l > 
and ( I26t . the result for bi is independent of the choice of 
N. Our numerical results for the two-loop Bethe logarithm 
of IS and 2S states read (results for 15 and 25 are quoted 
from El): 



(33c) 



b L (lS) - 


-81.4(3) , 


(34a) 


b L {2S) = 


-66.6(3) , 


(34b) 


b L {iS) = 


-63.5(6) , 


(34c) 


&l(45) - 


-61.8(8) , 


(34d) 


b L (5S) = 


-60.6(8) , 


(34e) 


b L (6S) = 


-59.8(8) . 


(34f) 



These results are displayed in Fig. [5] 

From here on, we restore in the following formulas the 
physical dimensions of all energies and frequencies and re- 
voke the scaling introduced in Eq. dl6l l. Primed quantities 
will no longer be used in the following sections of this work. 
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1 2 3 4 5 

principal quantum number n 



FIG. 5: (color online). Dependence of the two-loop Bethe logarithm 
&l (nS) on the principal quantum number n. The explicit numerical 
results [Eq. <34H are displayed together with a three-parameter fit 
of the form —57.4 — 13.7/n — 10.1/n 2 from which one may infer 
lim n _,oo bL(nS) — —57(2). Quantities plotted along the abscissa 
and the ordinate axis are (of course) dimensionless. 

IV. AMBIGUITY IN THE DEFINITION OF B eo 

Low Jl7ll was the first to point out that the definition of an 
atomic energy level becomes problematic at the order of a s 
[more specifically, a 2 (Za) 6 m], and that it becomes neces- 
sary at this level of accuracy to consider the contribution of 
nonresonant energy levels to the elastic scattering cross sec- 
tion. In lU9ll . it has been stressed that nonresonant effects are 
enhanced in differential as opposed to total cross section, lead- 
ing to corrections of order a 2 (Za) 4 m. Related issues have 
recently attracted some attention (see also the discussion in 
Sec.U, and there is even a connection to the two-loop cor- 



rections of order a 2 (Za) 6 m, as we discuss in the following. 
Namely, as pointed out in 11611 . the two-loop self-energy con- 
tains contributions which result from squared decay rates. 

For excited reference states, the nonrelativistic two-loop 
self-energy dl4t contains imaginary contributions which are 
generated by both the uj\- as well as the L02 -integrations. The 
imaginary part of the one-photon self-energy is generated by 
a pole contribution and leads to the decay rate which is the 
imaginary part of the self-energy. Consequently, real contri- 
butions to the two-photon self-energy which result as a prod- 
uct of two imaginary contributions are naturally referred to as 
squared decay rates. These are natural contributions to the 
two-loop self-energy shift in the order of a 2 (Za) 6 m and 
cannot be associated in a unique manner with one and only 
one atomic level. Roughly speaking, the problems in the in- 
terpretation originate from the fact that the Gell-Mann-Low- 
Sucher 1561 15711 formalism involves a priori asymptotic states 
with an infinite lifetime (vanishing decay rate). Furthermore, 
it has been mentioned 1 20] that problematic issues persist even 
if the concept of an atomic energy is generalized to a reso- 
nance with a finite width — i.e., even if the canonical concept 
of a pole of the resolvent on the second Riemann sheet 158115911 
is used for the definition of an atomic resonance. In general, 
the squared decay rates illustrate that we are reaching the limit 
of the proper definition of atomic energy levels in considering 
higher-order two-loop binding corrections. 

In fl6ll . the squared decay rates have been analyzed in some 
detail. There are four specific terms out of the nine in curly 
brackets in Eq. HAi which give rise to squared decay rates. 
We list these terms here with a special emphasis on th e hig her 
excited 35 state, following the notation introduced in 11611 . In 
the following formulas, the physical dimensions of all ener- 
gies and frequencies are restored [cf. Eq. 1161 1. and we have 
for the first term 71(35), which is the analog of Eq. (4) of ficii : 



Ti(35) = lim - 



2a 



s—>o+ V 3 7r m 



x ( 35 



H — iS — Ess + u>i H — Ess + lui + LU2 H — i 8 — Ess + ^2 



35 



(35a) 



For the analog of Eq. (8) of CI, 

we have 



T 2 (35) = lim - 



2a 



S— >o+ V 3 7r m 



x ( 35 



■P 



1 



H — i 8 — Ej,s +wi V H — £"35 / H — 18 — E 3S + u 2 



1 



P 



35 



(35b) 
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and we also have [see Eq. (15) of lll6ll l 



T 3 (35) = lim . 

5^0+ \ 3 7r m 1 



x ( 35 



H — iS — E%s + ^i 
The last relevant term is [see Eq. (17) of fl^l l 



Jo 



35 ) 35 



i? - i (5 - £3,5 + w 2 



35 



7^(35) = lim 



2a 



x ( 35 



dwi Wl / doj 2 ^2 

5^0+ V 37T771 2 / ./ J 



H — id — E3S +1^1 H — i S — £35 + 0J2 



35 



(35c) 



(35d) 



Here, 7? is the Schrodinger Hamiltonian. We now proceed 
to analyze the squared decay rates generated by the terms % 
(i = 1, ... ,4) in some detail. It should be reemphasized 
here that the main contributions to the energy shift gener- 
ated by the % have already been analyzed in Sec. [HI] How- 
ever, the prescriptions ( 1271 and d28i lead to a complete ne- 
glect of the (squared) imaginary contributions. Consequently, 
we here "pick up" only the terms of the "squared-decay" 
type — i.e., the terms generated by the infinitesimal half-circles 
around the poles at u)\ = E^s — E<ip and cj 2 = E$s — E-2p. 
For the evaluation of these squared pole terms, specification of 
the infinitesimal imaginary part — i 6 is required in order to fix 
the sign of the pole contribution. This procedure of extracting 
squared ima gina ry parts leads to the terms Cj (i = 1, . . . , 4), 
respectively l60ll . 

We now proceed to analyze the squared decay rates gen- 
erated by the terms % (i = 1, ... ,4) in some detail. The 
term T\ is due to the diagram with crossed loops in Fig.^A). 
For the contribution Ci(35) generated by the poles at u>\ = 

E^s — E2P and cj 2 = E 3 s — E2P in T\ (2P), we obtain 



d(3S) = a?^(E 3S -E 2P ) 2 \(2P\p\3S}\ 2 



x ( 2P 



H + E3S — 2 E2P 



2P 



trix element M 1 reads 



Mi 



1 



2P 



1 



■P 



H + E3S — 2 E2P 
and we have for the well-known dipole matrix element 



2P ) = 0.697, 
(37) 



2P 



35 



9 q3 



2^3 



(Zaf 



(38) 



Note that the contribution C\ lacks the factors 7r in the denom- 
inator which are characteristic of other two-loop corrections: 
these are compensated by additional factors of ir in the numer- 
ator that characterize the pole contributions. 

The rainbow diagram in Fig.[nB) with the second loop in- 
side the first does not create squared imaginary contributions. 
From the irreducible part of the loop-after-loop diagram in 
Fig-EC) (we exclude the reference state in the intermediate 
electron propagator), the term is obtained. Again, pick- 
ing up only those terms which are generated by the infinites- 
imal half-circles around the poles at lo\ = £35 ~ E2P and 
^2 = E3S — E2P, we obtain the contribution C2(35) involv- 
ing squared decay rates: 



C 2 (35) = a 2 



9m 4 



(E 3 s~E 2 pf |(2P|p|35)| 



x ( 2P 



1 



H — E 3 s 



2P 



3 3 5 ! 



■a 2 (ZafmM 1 



(36) 



where we define \2P) to be the state with magnetic quan- 
tum number (angular momentum projection) m = 0. This 
explains the additional factor of 3 in comparison to Eq. (5) 
of d. The factor originates from the summation over mag- 
netic quantum numbers of the \2P) state, and we reemphasize 
that we understand by \2P) only the state with magnetic quan- 
tum number (angular momentum projection) m = 0. Thema- 



2 5 
335* 



a 2 (Zaf m 



where the matrix element M2 reads 



Mi 



2P 



1 



H - E 3S 



(39) 



2P 



0.490. 



(40) 

The prime in the reduced Green function indicates that the 35 
state is excluded from the sum over intermediate states, and 
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TABLE I: Sample values of the g function, defined in Eq. <25> . for the nS states with 
n — 1, ... ,6. Multiplication by a factor of r? approximately accounts for the n-dependence, 
in agreement with the n~ 3 -type scaling of the two-loop correction as implied by Eq. 0. 





5is 


8g2S 


27 ff 3S 




125 g 3S 


216 g 6S 





0.000 00 


0.000 00 


0.0000 


0.0000 


0.0000 


0.000 


5 


-10.28160 


-10.367 94 


-10.4501 


-10.490 8 


-10.522 6 


-10.5460 


20 


-16.560 34 


-16.415 97 


-16.393 4 


-16.385 1 


-16.3861 


-16.3867 


80 


-22.71402 


-22.439 66 


-22.372 


-22.345 5 


-22.332 


-22.326 3 


180 


-26.232 35 


-25.923 09 


-25.848 


-25.813 6 


-25.798 1 


-25.789 5 



it should not be confused with the notation used in Sec. |in| 
where the prime was used to denote scaled dimensionless in- 
stead of dimensionful quantities. 

From the derivative term (reducible part of the loop-after- 
loop diagram), we obtain 



C 3 (2P) 



3m 
2 17 3 3 

K19 



I (E 3S -E 2P ) |(2P|p|35)r 



a 2 (Za) 6 m . 



(41) 



In order to derive the imaginary parts, one should remember 
that the squared propagator originates from a differentiation of 
a single propagator with respect to the energy. An integration 
by parts is helpful. 

The last contribution of the "squared-decay" type — it orig- 
inates from the "seagull term" characteristic of NRQED — is 
T 4 . The corresponding C-term is 



C* 4 (3S*) 



-a 2 — 
3m 3 



3 (E 3S - E 2P f |(2P|p|3S}| 2 



3 2 5 8 



a 2 (Za) e m . 



(42) 



Adding all contributions, we obtain a shift of 



4 

E 

i=l 



a 



(*,-(2) 



a\ 2 (Za) e m 



3 3 



(-0.00151) (43) 



for the 35 level. The numerical value is tiny; for Z = 1 the 
shift amounts to only 



5 2 v(3S) = -0.00565 Hz. 



(44) 



For the corresponding ambiguous contributions to the i?60 _ 
coefficient [see Eqs. and J43H . we use the notation 



S 2 B 6Q {3S) = -0.00151 



(45) 



For the 4S state, we have to take into account the decays into 
the 2P and 3P states. For example, the contribution Ci(4S) 
reads 



d(4S) = a 2 



9m 4 



(E 4S -E 2P ) 2 |(2P|p|4S)| 2 (2P 



1 



(E 4S -E 3P ) 2 |(3P|p|4S)| 2 (3P 



H + E 4 s — 2 E 2 p 
1 



2P 



P 



+a 



H + E4S — 2 E 3 p 
{E 4S - E 2P ) {E 4S - E 3P ) Re ( (4S\p>\2P) Up 



P 



3P 



9m 4 
a\ 2 {Zoifm 



1 



H + E 4 s — E 3 p — E 2 p 



3P) (3P|y'|4S) 



4 :i 



(0.00108) 



(46) 
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TABLE II: Squared decay rates are extracted as the squared bound-state pole terms from the terms 71 — To, in Eqs. I35al — J35dl . Explicit 
formulas (35 state) for the terms d (i = 1, . . . , 4) are given in Eqs. 1361 — 1421 . All contributions d scale as Z 6 , whereas the decay rates 
r given in the eighth column scale as Z 4 . Numerical values are given for Z = 1. The decay rates may be derived in the standard way 
[see the derivation of Fermi's golden rule as given in Eqs. (2. 103) — (2. 118) of loHl l. We only indicate approximate values for T, without 
relativistic corrections. For the 2P 1 / 2 states, a detailed calculation leads to T(2P 1 / 2 ) — 0.9970942 Z 4 MHz l54ll . For any given state, 
the squared decay rates 8 2 v are about seven to eight orders of magnitude smaller than the width F. All states listed in the table may decay 
via the El mode, wherefore the decay rates as well as the ambiguities 8 2 u are formally of the same order-of-magnitude [a (Za) 4 m and 
a 2 (Za) 6 m, respectively]. However, the numerical coefficients differ by two orders of magnitude; S states typically have a much longer 



lifetime. 


state 


Ci 




c. 


c 4 




8 2 B 60 


r 


T 




2P 


1.42208 Hz 


2.06790 Hz 


-1.00843 Hz 


-4.84593 Hz 


-2.36438 Hz 


-0.18789 


99.76 MHz 


0.16 x 10~ 8 


s 


3P 


0.50353 Hz 


0.06037 Hz 


-0.12787 Hz 


-2.00952 Hz 


-1.57349 Hz 


-0.42202 


30.21 MHz 


0.53 x 10" 8 


s 


3S 


0.00210 Hz 


0.00148 Hz 


-0.00018 Hz 


-0.00565 Hz 


-0.00564 Hz 


-0.00151 


1.00 MHz 


15.83 x 10~ 8 


s 


4.5' 


0.00170 Hz 


-0.00100 Hz 


-0.00015 Hz 


-0.01019 Hz 


-0.00964 Hz 


-0.00613 


0.70 MHz 


22.65 x 10~ 8 


s 



The sum of Ci 4 for the 45* state of hydrogenlike systems 
with (low) nuclear charge Z is 

ta(4S)=Q 2 <^ (-0.00613). (47) 

For atomic hydrogen (Z = 1), this correction evaluates to 

5 2 v(AS) = -0.00964 Hz. (48) 

This is numerically larger than the corresponding effect for 35* 
[see Eq. M l. We have 

S 2 B 6Q (45) = -0.00613. (49) 

Although self-energy corrections canonically scale as n~ 3 
[see Eq. Qj], the coefficient in this case grows so rapidly with 
n that the correction is enhanced for 45 in comparison to 35. 
Further detailed information can be found in Table [H] We 
also take the opportunity to clarify that numerical values for 
squared decay as given in fToll (for 2P and 3P states) should 
be understood as given in inverse seconds rather than Hz (see 
also the discussion near the beginning of Sec.lIIH. 

V. FURTHER CONTRIBUTIONS TO B m 

The coefficient B^q can be represented as the sum 

B w = b L + b M + b F + b H + by P . (50) 

The two-loop Bethe logarithm b^ comes from the region 
where both photon momenta are small and has been the sub- 
ject of this work. &m stems from an integration region where 
one momentum is large ~ m, and the second momentum is 
small. This contribution is given by a Dirac 5 correction to 
the Bethe logarithm [see also Eq. J2 It and Ref. 14211 1. It has 



already been derived in fl2il but not included in the theoretical 
predictions for the Lamb shift: 



b M = —N(nS). (51) 



As has already been mentioned in l55ll . the contributions bp 
and bff originate from a region where both photon momenta 
are large ~ m, and the electron momentum is small and large 
respectively. Finally, feyp is a contribution from diagrams that 
involve a closed fermion loop. None of these effects have 
been calculated as yet. On the basis of our experience with 
the one- and two-loop calculations we estimate the magnitude 
of these uncalculated terms to be of the order of 15%. For 
higher excited states (35, . . . , 6S), the 15 % uncertainty due 
to unknown contributions is larger than the ambiguity 8 2 Bqq 
listed in Table |TT| 

Concerning logarithmic two-loop vacuum-polarization ef- 
fects fl2ll . we mention that the contribution of the two-loop 
self-energy diagrams to Bqi for the IS state reads 49.8, 
whereas the diagrams that involve a closed fermion loop 
amount to 0.6. Concerning the one-loop higher-order binding 
correction Aqq(1S) (analog of Bqq) it is helpful to consider 
that the result for 15 is -30.92415(1) (see Refs. SEl); this 
is the sum of a contribution due to low-energy virtual photons 
of —27.3 fl2l Eq. (5.116)], and a relatively small high-energy 
term of about —3.7 ILa Eq. (6.102)]. In estimating these con- 
tributions, we follow 15511 . 

This leads to the following overall result for the Bqq co- 
efficients, where the first two results (15 and 25) are quoted 
from l55ll . and the latter results are obtained within the current 
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investigation: 



Beo{lS) = 


-61.6(3) ±15%, 


(52a) 


B 60 (2S) = 


-53.2(3) ±15%, 


(52b) 




-51.9(6) ± 15% , 


(52c) 


B 60 ^S) = 


-51.0(8) ±15%, 


(52d) 


B 60 (5S) = 


-50.3(8) ±15%, 


(52e) 


B 60 (6S) = 


-49.8(8) ± 15% . 


(52f) 



The values are in numerical agreement with those used in 
latest adjustment of the fundamental physical constants 0]; 
these are based on an extrapolation of the results obtained for 
n = 1,2 l55Tl to higher n, using a functional form a+b/n, 
with an extra uncertainty added in order to account for the 
somewhat incomplete form of the functional form used in the 
extrapolation. We here confirm the validity of the approach 
taken in |8] by our explicit numerical calculation. 



VI. CONCLUSIONS 

The calculation of binding two-loop self-energy corrections 
has received considerable attention within the last decade l26l 
I27I Eijl l62ll . As outlined in Sec. U there is an intuitive 
physical reason why a reliable understanding of the two- 
loop energy shift requires the calculation of all logarithmic 
as well as nonlogarithmic corrections through the order of 
a 2 (Za) 6 m. It is the order of a 2 (Za) 6 m which is the "nat- 
ural" order-of-magnitude for the two-loop self-energy effect 
from the point of view of nonrelativistic quantum electrody- 
namics (NRQED); i.e., low-energy virtual photons begin to 
contribute at this order only, whereas effects of lower order 
[a 2 (Za) 4 m and a 2 {Za) 5 m] are mediated exclusively by 
high-energy virtual quanta. 

In Sec. [n] we recall known lower-order coefficients for S 
states, as well as logarithmic corrections. The formulation 
of the problem within NRQED l5(ill and the actual numeri- 
cal evaluation of the two-loop Bethe logarithms for higher 
excited S states is discussed in Sec. [Hi] Numerical results 
for 6^ are given in Eq. J34b . As shown in Fig. [5J the de- 
pendence of these results on the principal quantum number 
follows a pattern recently observed quite universally for bind- 
ing corrections to radiative bound-state energy shifts I42ll63ll . 
This permits an extrapolation of the results to higher princi- 
pal quantum numbers, which is useful for the determination 
of fundamental constants 1 8 ] . 

There is a certain ambiguity in the definition of the two- 
loop nonlogarithmic coefficient Bqq due to squared decay 
rates (Sec. IIVI : this aspect has previously been considered 
in fl6ll for P states. Here, the treatment of the squared decay 
rates is generalized to excited S states. The ambiguity, while 
formally of the order of a 2 (Za) 6 m, is shown to be barely 
significant for S states (see TablellH. due to small prefactors. 



Numerical estimates of the total Bqq -coefficient for excited 
nS states (n = 1, . . . , 6) are given in Eq. d52i . These results 
improve our theoretical knowledge of the hydrogen spectrum. 
On the occasion, we would also like to mention ongoing ef- 
forts regarding the calculation of binding three-loop correc- 
tions of order a 3 (Za) 5 @|. At Z = 1, these binding three- 
loop corrections are of the same order-of-magnitude (a 8 ) as 
the two-loop Bethe logarithms discussed here. There is con- 
siderable hope that in the near future, our possibilities for 
a self-consistent interpretation of high-precision laser spec- 
troscopic experiments may be enhanced significantly via a 
combination of ongoing experiments at Paul Scherrer Institute 
(PSI), Laboratoire Kastler-Brossel and Max-Planck-Institute 
for Quantum Optics, whose purpose is a much improved 
Lamb-shift measurement (15-25'- and 15'-3S'-transitions 
combined with an improved knowledge of the proton charge 
radius as derived from the PSI muonic hydrogen experiment). 
The comparison of numerous transitions in hydrogenlike sys- 
tems with theory may also help in this direction as it allows 
for an evaluation of the proton charge radius using an overde- 
termined system of equations, provided that theoretical Lamb- 
shift values are used as input data for the systems of equations 
rather than variables for which the systems should be solved 
[see, e.g., Eqs. (2) and (3) of P i. 

Finally, we would like comment on the relation of the an- 
alytic approach (Za expansion) pursued here and numeri- 
cal calculations of the self-energy at low nuclear charge Z, 
which avoid the Za expansion and which have been car- 
ried out on the one-loop level for high nuclear charge num- 
bers Z 16611 . with recent extensions to the numerically more 
problematic regime of low Z [5| . One might note that tradi- 
tionally, the most accurate Lamb-shift values at low Z have 
been obtained via a combination of analytic and numerical 
techniques — i.e., by using both numerical data obtained for 
high Z and known analytic coefficients from the Za expan- 
sion l67ll . (This is one of the main motivations for pursuing 
both numerical and analytic calculations of the two-loop self 
energy, in addition to the obvious requirement for an addi- 
tional cross-check of the two distinct approaches.) The gen- 
eral paradigm is the extrapolation of the self-energy remain- 
der function obtained from high-Z numerical data after the 
subtraction of known analytic terms; in many cases this ex- 
trapolation leads to more accurate predictions for the remain- 
der at low Z than the simple truncation of the Za expansion 
alone. Various algo rithms have been developed for this pur- 
pose (see, e.g., l67ll68Tl ). Indeed, the combination of analytic 
and numerical techniques has recently proven to be useful 
in the context of binding corrections to the one-loop bound- 
electron g factor |69], although direct numerical evaluations 
at Z = 1 had been available before |7|]. Still, it was possi- 
ble to improve the theoretical predictions for the g factor self- 
energy remainder function at low Z by an order of magnitude 
via a combination of the analytic and numerical approaches, 
in addition to the fact that an important cross-check of the an- 
alytic and the numerical approaches versus each other became 
feasible. 
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